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Abstract. We present a class of symplectic matrices which transform by similar- 
ity given 2n x 2n -dimensional matrix into Bunse-Gerstner form. If the given matrix is 
skew-Hamiltonian, the transformation gives a solution of an antisymmetric Riccati matrix 
equation, resulting from optimal control of linear systems. 
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1. Introduction. 



The paper is concerned with the class of 2n x 2n dimensional matrices: 

S- 



Su 
S21 



S12 
S22 



XI) 



where Sij G {1, 2}) are arbitrary real n x n matrices. Let the matrix J be defined by: 



J 



In 





A matrix S is Hamiltonian if J T SJ = —S T . The Hamiltonian matrices come out from 
optimal control of linear systems [1],[2],[4],[5],[7]-[11],[17]. A matrix S is skew-Hamiltonian 
if J T SJ = S T . For example, if S is a Hamiltonian matrix, then the matrix S 2 is skew- 
Hamiltonian. 

We represent some further definitions and basic results. A matrix U is symplectic if 
U T JU = J. The symplectic similar transformations keep the Hamiltonian structure of 
matrices, i.e. if S is a Hamiltonian matrix and U is a symplectic matrix, then U~ 1 SU is 
also a Hamiltonian matrix. Analogously, it is easy to verify that the symplectic similar 
transformations also keep the structure of the skew-Hamiltonian matrices. The Hamilto- 
nian matrices form Lie algebra and the symplectic matrices form the corresponding Lie 
group. 

Bunse-Gerstner proved [6] that by application of an orthogonal symplectic transfor- 
mation one can reduce an arbitrary real matrix S (1.1) in a form 



SI Q/ 
11 °12 

J 21 °22 



(1.2) 
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such that its (2,1) block S' 21 is upper-triangular and (1,1) block S'^ is upper-Hessenberg. 
For the transformation matrix a composition of finite number of Givens-like [16] and 
Hauseholder-like [16] symplectic orthogonal transformations is used. 

Van Loan proved [15] that if S is a skew-Hamiltonian (squared Hamiltonian) matrix, 
then S is a symplectic-orthogonal similar to a matrix of the form: 

\S' n s> 



where S' n is a Hessenberg matrix and S' 12 is a skew-symmetric matrix. For the transfor- 
mation matrix the same composition of finite number of Givens-like and Hauseholder-like 
symplectic orthogonal transformations is used. Besides Van-Loan's method [15], [3], [6], for 
this purpose, the methods from [4] or [13] may be applied. 

In [14] are applied Arnoldi-like matrices for similarity transformation of 2n x 2n - 
dimensional matrix S into a form (1.2), such that the submatrix S 21 = 0, under the 
condition that a minimal annulator of the first column of the transformation matrix, with 
respect to S, is of degree n. Although for general matrices S the last condition is connected 
with solving a nonlinear eigenvalue-eigenvector problem, for the skew-Hamiltonian matrices 
S these transformations result in S' 2 i = also, as do in Van Loan's work [15], because the 
minimal annulator of skew-Hamiltonian matrices is of degree n. 

In this paper we present a new class of symplectic matrices which transform 2n x 2n 
-dimensional matrices into Bunse-Gerstner form (1.2), instead of the Givens-like, Haush- 
older -like or Arnoldi-like symplectic transformation matrices. We use the matrices of the 
form: 

(1.4) 

where Q is an orthogonal matrix and Y is a symmetric matrix. This form arises naturally 
in solving of the antisymmetric Riccati matrix equation (ARME) [13]: 

-YS 12 Y + S 22 Y - YSu + S 21 = 

where S is a skew-Hamiltonian matrix and Y is an unknown symmetric matrix, because 
the matrix Y in (1.4) is already a solution of ARME, in case the matrix (1.4) brings S in 
block-triangular form. To compare our result with [15], if a solution of ARME is required 
after the Van Loan's transformation [15] is performed, one has to compute the matrix 
U 2 iU{[ , where U\\ and U 2 \ are blocks of the transformation matrix. If the matrix J7n is 
singular, it can not be done. 

Other references are mentioned in the text. 

2. Main results. 

Let be given the following 2n x 2n block matrix 

(2.1) 



(1.3) 



U 



Q 
YQ Q 
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where SV, G {1,2}) are arbitrary n x n matrices. In this section we present a finite 
algorithm for similar transformation by a matrix U of the matrix S, i.e. 



u~ 1 su = 



Of Qf 



Of 
D 21 



qf 



= S' 



where S'n is an upper Hessenberg matrix and S' 21 is an upper triangular matrix. 
Let us consider the following set of regular 2n x 2n matrices: 



G = {U 



L 
YL L~ T 



Y}. 



It is easy to verify that the set G with the matrix multiplication is a group, a subgroup of 
the group of the symplectic matrices. 

Applying the similar transformation U~ 1 SU, we obtain 



u~ 1 su = 



L- 1 
-L T Y 





L T 



S21 



S\2 
S22 



L 
YL L~ T 



L-^Sn + S^L L- 1 S 12 L~ T 
L T (S 2 i + S 22 Y - YSu - YS 12 Y)L L T (S 22 - YS 12 )L~ T 



Qf Of 

°11 D 12 

V 



S'. 



. S'21 S'22 . 

The matrix S will be transformed n — 1 times with respect to the matrices 



YiLi 







Y T — Y 



(1 < i < n- 1). 



(2.2) 



(2.3) 



By the first transformation, the first column of the new matrix S'21 will become upper 
triangular and the first column of the new matrix Sn will become Hessenberg. Continuing 
this process, by the i-th transformation, the i-th column of the new matrix S'21 will become 
upper triangular such that the previous % — 1 columns of that matrix will remain upper 
triangular, and the i-th column of the new matrix S\\ will become Hessenberg, such that 
previous % — 1 columns of that matrix will remain Hessenberg. 
From 



U = 
we obtain 



Li 







L 2 
Y 2 L 2 







L n -i 
Y n -\L n _i 







J n-1 





' L 







(2.4) 



L = L\L 2 

Y = Y X + L^ T Y 2 L^ 1 + L^L^YsL^Li 1 + ■ 



lj n-2 1 n-l^ n - 2 



L, 2 L, 1 



and hence Y is a symmetric matrix. Further we will show explicitly the structures of 
Li and Yi which are needed for defining the i-th step of the algorithm. For the sake of 
simplicity, the indices of the matrices Yi will be omitted. 
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Let us put 



Y = aw = a 





vv T 



Vl 



V = 



Vi+1 



(2.5) 



L' 



h 
L 



L = 



1 

W i+2 1 

w i+3 1 



w n 



lJ 



w 







- ' 


Wi+1 







W i+2 






- w n . 







where a is a scalar, v and w are vector-columns. If the i-th column of S21 is already upper 
triangular, we put Y = 0. The matrix Y is symmetric. 
Let us denote 



Sn = [n, . . . ,r n ] = 



rn ••• n, 



, = [ti, . . . ,t n ] = 



t 



11 



t 



-t n l 



In 



tnn - 



The i-th column of the matrix S' 2 \ is upper triangular, if and only if 



U — aw Ti = 



i.e. the last n — i elements of this vector should be equal to zero. If 

n 

j=i+l 

the solution of this system of a and v yields to 

Vj = tji , (i + 1 < j < n) and a 



(2.6) 



The vector w which takes part in the construction of the matrix L' , will be determined 
from the condition that the i-th column of the matrix L~ l SnL be Hessenberg. It will be 
done as follows. If that column is already Hessenberg, we put L = I. If it is not Hessenberg, 
then among the elements r^+i^, ri+2,1, ■ ■ ■ > fni we choose that one of maximal module. Let 



4 



Tki be such an element. For such k we do a similar transformation of the matrix Sn which 
consists of interchanging the k-th with the (i + l)-th row and also the k-th with the (i + 1)- 
th column of Sn. Thus the maximal module element will come on position {i + 1, i). For 
the sake of simplicity, the obtained matrix will also be denoted by Su, and the matrix L' 
which will be obtained further should be multiplied from left with the permutation matrix 
in order to obtain the matrix L. The permutation matrix on S belongs to the group G. 
The i-th column of L^SuL and hence of S'n is Hessenberg if and only if: 



1 





• 


■ o- 






-W i+2 


1 


• 


• 











1 • 


• 




r i+2,i 
S 

- fni 


-w n 





■ 


• 1. 







Hence the solution for the vector w is given by 



W i+2 = 



W n = 



r i+l,i 



r i+l,i 



We should prove that the first i — 1 columns of the new obtained matrix £21 remain 
upper triangular after the i-th transformation, assuming that the first i — 1 columns of the 
matrix S21 are upper triangular. Indeed, 

S'21 — L T S2iL + aL T S22VV 1 L — aL T vv T LL~ X SuL — a 2 v T S12V L T vv T L . 

Since the first i elements of the vector row v T L are equal to zero, the first i columns of the 
second and the fourth matrices in the previous sum, are zeros. The first i — 1 columns of 
the matrix L~ 1 SnL are Hessenberg according to the inductive assumption, and the i-th 
column has been done Hessenberg using the vector w. Thus the first i — 1 elements of the 
following vector-row 

v T LL~ 1 S 11 L 

are zeros. Consequently, the first i — 1 columns of the matrix S' 2 i will be upper triangular. 
The i-th column of that matrix is already upper triangular, using a and v. 

We should prove that the first i columns of the matrix S'n are Hessenberg. It follows 
from the following representation 

S' n = L^SuL + aL- 1 S 12 vv T L 

and from the fact that the first i columns of the matrix L~ 1 SnL are Hessenberg. 
According to the following forms of the matrices Li and Y t 
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it follows that the first row and the first column of the matrix Y are zero. 
We have proved the following theorem: 

Theorem 2.1. Let be given a matrix S by (2.1). Applying the similar transformations on 
S by the matrices Ui, % = 1, . . . , n — 1 (2.3), under the conditions (2.6), we get the global 
transformation matrix U (2.4) with matrix Y in which the first row and column are equal 
to zero, and the transformed matrix S' be such that S'n is an upper Hessenberg matrix 
and S' 21 is an triangular matrix, m 

Remark. One feature of the presented algorithm is that it contains only linear algebraic 
operations: adding, subtracting, multiplying and dividing. Thus the theorem is true also 
if the elements of the matrix S are elements of larger set. For example it can be the field 
of complex numbers or the ring of the analytical functions of more variables. 

However, if the elements of the matrix S are real numbers, from numerical viewpoint 
it is more convenient if the matrix L is orthogonal. Therefore, we define the following 
algorithm. 

For arbitrary given matrix S, we shall construct the matrix: 



U 



Q 
YQ Q 



Y 



(2.7) 



such that 



Q T 
-Q T Y Q T 



$11 S12 

$21 S22 



Q 
YQ Q 



0/ 0/ 

Of 0/ 
D 21 °22 



where S'n is an upper Hessenberg matrix and S f 21 is an upper triangular matrix, and the 
first row and column of Y are zeros. 

It verifies that the set of the matrices U with respect to the matrix multiplication is 
a group. We apply n — 1 such similar transformations to S determined by 



Ui 



Q l 
YiQi Qi 



Y? 



Yi 



(1 < % < n- 1). 



(2.8) 



By the first transformation, the first column of the obtained matrix S 2 i will become up- 
per triangular, and the first column of the obtained matrix Sn will become Hessenberg. 
Generally, applying the i-th transformation, the i-th column of the obtained matrix £21 
will become upper triangular such that the first i — 1 columns will remain upper triangular 
and the i-th column of the obtained matrix Sn will become Hessenberg such that the first 
i — 1 columns of that matrix will remain Hessenberg. 
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The structures of the matrices Yi are the same as in the proof of Theorem 2.1, and 
the matrices Qi without indices, for simplicity, are given by 



Q = I n — 2ww T 







I - 2ww 



w 



Wi 




' o ■ 


Wi 







Wi+1 




w i+ i 


- w n . 







, w 



W r 



T -T - 

w w = w w 



1, 



(2.9) 



where w is a vector-column whose first % elements are zeros. The matrix Q is Householder's 
[16] and orthogonal. 

The vector w which takes part in the construction of the orthogonal matrix Q, will 
be determined such that the i-th column of the matrix Q T SuQ to be Hessenberg, i.e. 





~ r i+l,i~ 






2ww T ) 


r i+2,i 













- fni - 




.0. 



(2.10) 



The system (2.10) together with the orthogonality condition (2.9) has the following solution 
for w: 

w i+ i = (3{r i+1)i + s) 



w i+ 2 = f3r i+2 ,i 



w n = f3r ni 

where 

(3 = [2(« 2 + r l+1 ,s)]-^ , 8 = [rl hi + + s + r^. 

We should prove that the first i — 1 columns of the new obtained matrix S21 remain upper 
triangular after the i-th transformation, assuming that the first i — 1 columns of the matrix 
S21 are upper triangular. We have 

S'21 = Q T S 2 iQ + aQ T S 22 vv T Q - aQ T vv T QQ T SuQ - a 2 v T S 12 vQ T vv T Q. 

Since the first i elements of the vector-row v T Q are zeros, it follows that the first i columns 
of the second and the fourth matrices in the previous sum are zeros. The first i — 1 columns 
of the matrix Q T SnQ are Hessenberg according to the inductive assumption, and the i-th 
column was made Hessenberg using the vector w. Thus the first i — 1 elements of the 
following vector-row 

v T QQ T SnQ 
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are zeros. Consequently, the first % — \ columns of the matrix S' 2l will be also upper 
triangular. The i-th column of that matrix was done upper triangular using a and v. 

The first row and column of the matrix Y are zero, from the same reasons as in the 
proof of Theorem 2.1. 

We have proved the following theorem. 

Theorem 2.2. Let be given a matrix S with (2.1). Applying the similar transformations 
on S by the matrices Ui, % = 1, . . . , n — 1 (2.8), under the conditions of Gniteness of f3 in 
each step, we get the global transformation matrix U (2.7) with matrix Y in which the 
first row and column are equal to zero, and the transformed matrix S' be such that S'n is 
upper Hessenberg matrix and S' 21 is upper triangular matrix. ■ 

As a consequence from these results follow two corollaries for the Hamiltonian and 
skew-Hamiltonian matrices. 

Corollary 2.3. Let the matrix S be a Hamiltonian matrix. Applying n — 1 similar 
transformations (2.3) of the matrix S under the conditions (2.6), we get a matrix S' , such 
that the left lower block of dimension n x n is a diagonal matrix and the upper left block 
of dimension n x n is of Hessenberg form. 

Proof. It follows from (2.2) that the similar matrix S' = U~ 1 SU also has the Hamil- 
tonian form, if U is any matrix of the form given in Theorems 2.1 and 2.2. The rest of the 
proof is obvious. ■ 

The second corollary is about the ARME. 

Corollary 2.4. Let S be a skew-Hamiltonian matrix. Applying n — 1 similar transfor- 
mations (2.3), under the condition (2.6), the matrix S can be reduced in a form such that 
lower left block of dimension n x n is zero, the upper left block of dimension n x n has a 
Hessenberg form. If the final similar transformation is given by the matrix U , i.e. 

Q 0" 
YQ Q\ ' 

then the matrix Y is a symmetric solution of the ARME 

-YS 12 Y + S 22 Y - YS n + S 21 = 0, 

such that its first row and column are zeros. 

Proof. Using the similar transformation S' = U^SU, where U is any of the matrices 
given in Theorems 2.1 and 2.2, it is obtained skew-Hamiltonian matrix 5". The rest of the 
proof is obvious. ■ 

Remark. Breaking down the algorithm happens when there is no solution of ARME 
or there is no solution of ARME with first row and column equal to zero. 

Remark. If the algorithm breaks down or if it is requested a solution of ARME which 
is non-singular matrix, or the matrix with given first row and column, we could introduce a 
new unknown symmetric matrix X by Y = M+N~ T XN~ 1 , where M is a given symmetric 



U = 



L 

YL 







or U = 



8 



matrix and N is a given non-singular matrix. The new equation of X is also of the ARME 
type: 

—XS 12 X — XS n + S 22 X + S 2 ± = , 



where 



S' u =N- 1 (S 11 + S 12 M)N = S! 2 T 2 , 



S[ 2 = N-'LSmN-' 1 = -S 



/T 
'12 5 



S' 21 = N T (S 2 i + S 22 M - MSn - MS 12 M)N = -S: 



21 ) 



S' 22 = N T (S 22 - MS l2 )N~ T = S£ 

(compare with (2.2).). 

Example. In this example it is n = 6 and the matrix S is skew-Hamiltonian with 
randomly chosen numbers from (0,1) and symmetrically from (—1,0). The beginning 
matrix S is 



5-11 = 



0.001251 
0.350291 
0.710501 
0.147313 
0.008911 
0.166234 



0.563585 
0.895962 
0.513535 
0.165899 
0.377880 
0.663045 



0.193304 
0.822840 
0.303995 
0.988525 
0.531663 
0.450789 



0.808740 
0.746605 
0.014985 
0.445692 
0.571184 
0.352123 



0.585009 
0.174108 
0.091403 
0.119083 
0.601764 
0.057039 



0.479873 
0.858943 
0.364452 
0.004669 
0.607166 
0.607685 



, S22 = -Sn 



Sl2 — 



- 0.000000 


0.783319 


0.802606 


0.519883 


0.301950 





875973" 


-0.783319 


0.000000 


0.726676 


0.955901 


0.925718 





539354 


-0.802606 


-0.726676 


0.000000 


0.142338 


0.462081 





235328 


-0.519883 


-0.955901 


-0.142338 


0.000000 


0.862239 





209601 


-0.301950 


-0.925718 


-0.462081 


-0.862239 


0.000000 





779656 


.-0.875973 


-0.539354 


-0.235328 


-0.209601 


-0.779656 





000000. 


' 0.000000 


0.843654 


0.996796 


0.999695 


0.611499 





392438" 


-0.843654 


0.000000 


0.266213 


0.297281 


0.840144 





023743 


-0.996796 


-0.266213 


0.000000 


0.375866 


0.092624 





677206 


-0.999695 


-0.297281 


-0.375866 


0.000000 


0.056215 





008789 


-0.611499 


-0.840144 


-0.092624 


-0.056215 


0.000000 





918790 


.-0.392438 


-0.023743 


-0.677206 


-0.008789 


-0.918790 





000000. 



S21 — 



After the similar transformation (2.2), the matrix S becomes 



0.001252 
-0.822757 
0.000000 
0.000000 
-0.000000 
-0.000000 



2.394542 
0.705668 
4.483943 
-0.000000 
-0.000000 
-0.000000 



0.628127 
-0.064203 

1.443984 
-0.232568 

0.000000 
-0.000000 



-1.333712 
-0.380534 

3.651134 
-1.192501 
-2.882063 

0.000000 



1.170301 
0.078576 

-1.359573 
0.923134 
1.947789 

-0.112915 



-2.988671 
-2.116390 
-5.058558 
-0.061005 
-0.289124 
-0.049842 



S'22 = Sn 
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s 



12 



0.000000 
1.299940 
0.502376 
0.391241 
-0.293239 
-0.453059 



-1.299940 
-0.000000 
-1.315821 
0.067147 
-0.335201 
-0.359848 



-0.502376 
1.315821 
0.000000 

-0.003709 
0.814256 

-1.134890 



-0.391241 
-0.067147 

0.003709 
-0.000000 
-0.579926 

0.053570 



0.293239 

0.335201 
-0.814256 

0.579926 
-0.000000 

0.078375 



0.453059 ■ 
0.359848 
1.134890 
-0.053570 
-0.078375 
0.000000 . 



S 2 i = 

The symmetric solution of ARME with zero first row and column is 



Y 



-0.000000 
-0.186413 
-1.002085 
-0.065432 
-2.829218 

L -0.000000 -0.186413 -1.002085 -0.065432 -2.829218 2.524725 J 



-0.000000 0.000000 -0.000000 0.000000 0.000000 

0.000000 -2.976399 0.633676 -1.401061 -1.558562 

-0.000000 0.633676 -1.373721 -0.563505 0.769165 

0.000000 -1.401061 -0.563505 -0.599371 -1.051505 

0.000000 -1.558562 0.769165 -1.051505 1.478251 



Conclusions. We have presented a class of symplectic matrices which we have used 
in transformation of matrices S into Bunse-Gerstner form. Our numerical experiments 
with skew-Hamiltonian matrices S showed that for larger n (n > 10), the elements of the 
matrix S' 21 are not zeros, as theoretically they have to be. Further research could be to 
find the reasons for that behavior and to find an improvement of the presented algorithm 
for annihilating the (2,1) block of skew-Hamiltonian matrices. For example, one could try 
with multiple application of the algorithm. 
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